The ID t—J model with next-nearest neighbor hopping - breakdown of the Luttinger 

liquid? 
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We investigate the effect of a next-nearest neighbor hopping integral t 1 in the ID t—J model, 
using Lanczos diagonalization of finite chains. Even moderate values of t' have a dramatic effect 
on the dynamical correlation functions and Fermi surface topology, the high-energy holon bands 
become diffuse and overdamped, the band structure near the Fermi energy is dominated by t' . With 
increasing hole concentration the system underdoes a phase transition which changes the volume of 
the Fermi surface, but neither phase has a k,F compatible with a Luttinger liquid. 
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One dimensional (ID) systems of interacting Fermions 
have been the subject of much interest over many years. 
Based on the results of exact solutions for some models 
and approximate techniques, such as Bosonization 
[ 3 1 or the renormalization group || , the belief has formed 
that the universal fixed point model for systems with a 
repulsive interaction between the Fermions is the Lut- 
tinger liquid. Key features of this exotic quantum liquid 
are the existence of 'noninteracting' collective spin and 
charge like excitations, the so-called spinons and holons, 
which usually have different velocities, and the two Fermi 
points which obey the Luttinger sum-rule and govern the 
low energy physics. In the following, we want to present 
numerical evidence indicating that a very simple (and 
actually quite "physical") model, the ID t-J model with 
next-nearest neighbor hopping, may not fall into the class 
of the conventional Luttinger liquid. More precisely, the 
model reads: 



H 



tc) 



CT Cj±l,cr 



t c! Ci ±2)I 



Si 



Here Cj i(T = c 



(1 — rii^g) and Si denotes the spin oper- 
ator on site i. Throughout we keep the values t—l and 
J=0.2 (although the results do not depend on this spe- 
cial choice). 

In the "ordinary" t-J model the combination of hard-core 
constraint and one-dimensionality obviously does not al- 
low the sequence of spins along the chain to be changed 
by the motion of holes. The factorization of the wave 
function into a charge and spin part [pj then appears as a 
quite natural consequence. Introduction of next-nearest 
neighbor hopping changes this, in that hopping alone 
now can interchange spins. Obviously, this re- introduces 
some coupling between spin and charge degrees of free- 
dom, but the question is whether this coupling does not 
simply "renormalize to zero" in the actual ground state. 
We tried to address this question by studying various dy- 
namical correlation functions, computed numerically by 
Lanczos diagonalization of small clusters 0). It turned 
out that even for quite moderate values of t' the cou- 



pling of spinon and holon is large and in fact seems to 
be strong enough to induce a breakdown of the Luttinger 
Fermi surface. 




FIG. 1. Electron removal spectrum for a half- filled 20-site 
ring with t' = (left) and t' = -0.4 (right). The dotted fine 
indicates the "holon band" . 



To begin with, we consider the electron removal spec- 
trum at half-filling shown in Figure Introduction of 
even a rather small t' has a quite significant effect on the 
spectrum. For t'=0 the spectrum consists of a sequence 
of sharp peaks, which actually form a very systematic 
network of "spinon and holon bands" fffl]. The topmost 
band for k<n/2 traces out the spinon dispersion, and 
consequently its bandwidth scales strictly with J, the 
topmost band for k>n/2 follows the holon dispersion M. 



1 



In the spectrum for t'j^O some diffuse remnant of the 
"holon band" can be identified, but it no longer has the 
character of a well-defined excitation. The smaller peaks 
at higher excitation energies, which were sharp and had 
a well-defined dispersion for t'—O, disappear completely. 
Obviously the coupling of spin and charge degrees of free- 
dom due to t' induces a strong damping of the holon. 
Next, the uppermost band between k=0 and k=ir/2 has 
inverted dispersion, so that its top now is at fc=0. A more 
detailed study of this "quasiparticle band" for different t' 
(see Figure ^) shows that its dispersion is determined by 
t'. The width of this band, defined as W = e(0) - e(ir/2) 
can be fitted well by the expression W(t')~\.2 ■ (J+t'). 
Its dispersion on the other hand seems to change from 
nearest neighbor hopping, e(k) oc cos(fc) to next-nearest 
neighbor hopping, e(k) oc cos 2 (fc) for larger t'. This 
shows that the character of this band must be very dif- 
ferent now in that its 
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FIG. 2. "Quasiparticle band" in the electron removal spec- 
trum of a half- filled 20-site ring with different values oft'. The 
dotted line gives the curves e(fc) = W(t') cos(fc) for t' = 0, 0.2, 
and e(k) = W(t') cos 2 (fc) for t' = 0.4, 0.6 

dispersion is determined by a combination of the ex- 
change constant (which acts only on the spin degrees of 
freedom) and the next-nearest neighbor hopping integral 
(which acts on the charge degrees of freedom). Already 
this simple comparison thus demonstrates that the t'- 
term has a surprisingly strong effect. 
We proceed to the doped case and study the spin and 
charge-density correlation function, shown in Figure 
The density correlation function (DCF) shows a striking 
difference between the case t'=0, where it consists of a 
series of sharp dispersive peaks, and the case t'=— 0.4 
where the peaks are replaced by structureless continua. 
It is only at very low excitationen energies that there are 
sharp peaks also for nonvanishing t' ', the lowest of them 
appearing at fco.c=27r/8. Again, we see the strong damp- 
ing of the holon introduced by t' . On the other hand, 



the spin correlation function (SCF) does show some well 
defined peaks. A major difference is the very small en- 
ergy scale of the SCF for t'=— 0.4. After rescaling ex- 
citation energies, however, there is a certain similarity 
with the SCF for t'=0. One can still recognize a sc- 
ries of faint peaks which trace out the "spinon arc" 
and take their minimal excitation energy at 2kp =7tt/8. 
These peaks, however, correspond in fact to relatively 
high excited states, and the lowest peak, which also has 
a much higher intensity, now occurs at ko^—n/8. As 
would be the case in a conventional Luttinger liquid 
we thus have fco,s=&o,c/2, which interpretation however 
would force us to choose 2fcp=7r/8, corresponding to a 
Fermion density of 1/8, i.e. the density of holes in the 
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FIG. 3. Dynamical spin and charge correlation function 
for a 16-site ring with two holes with t'=0 (dotted line) and 
t'=— 0.4 (full line) Excitation energies for the SCF at t'=0 are 
multiplied by 0.12. 

half-filled band. To further address this surprising re- 
sult, we proceed to the complete single-particle spectral 
function, shown in Figure ^ for different hole number. 
We begin with the case of two holes (left hand panel) 
and focus on energies around Ep. There, the "quasi- 
particle band" persists with nearly unchanged spectral 
weight and dispersion, and the chemical potential essen- 
tially cuts into this band to form a hole pocket at fc=0. 
While the qp-peak at fc=0 crosses completely to the in- 
verse photoemission side, there is also some low intensity 
IPES weight at fc=7r/8, which probably corresponds to 
the standard "smearing" of the Fermi surface due to in- 
teractions. Doping two holes into the system shifts one 
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/c-point (namely fe=Q) through the chemical potential; 
this implies that the doped holes act as spin- 1/2 parti- 
cles, leading to the "effective" Fermion density 1/8, 
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FIG. 4. Full single particle spectral function for a 16-site 
ring. The chemical potential, defined as being halfway be- 
tween the excitation energies of first ionization and first affin- 
ity state, is the zero of energy, t = —0.4. 

in complete agreement with the hp inferred from the cor- 
relation functions. On the other hand, there are also 
deviations from pure rigid-band behaviour: starting at 
approximately 37r/8 there appears a broad high-energy 
band at an excitation energy of approximately 0.5t in 
the inverse photoemission spectrum. This "band" has 
virtually no dispersion, but is spectral weight increases 
continuously with k. It never even approaches the Fermi 
energy and thus is obviously unrelated to any "Fermi 
surface physics" (this is also confirmed by the ground 
state momentum distribution, which decreases slowly 
and completely continuously for k>2ir/8). It should be 
noted that the appearance of considerable electron ad- 
dition weight in the outer part of the Brillouin zone is 
an essentially trivial property of any Hamiltonian with 
a dominant nearest-neighbor hopping term: it is simply 
a necessary condition for having negative kinetic energy 
. In the 2D model this high energy inverse photoemis- 
sion weight in the outer half of the Brillouin zone can 
be explained by "charged magnons" , i.e. spin excitations 
which previously formed the dressing could of the an- 
nihilated hole ||. It is quite plausible that the diffuse 
high energy band in the present case is of a similar ori- 
gin. While the spectral function for two holes thus shows 



rather clear rigid-band behaviour (at energies close to 
Ep) the situation changes completely for 4 holes (see the 
right panel of Figure ||) . To begin with, unlike the two- 
hole case where A(k,u>) had extended incoherent con- 
tinua, the spectral weight is now concentrated in very 
few sharp peaks (for graphical reasons we have used dif- 
ferent Lorentzian broadenings and y-axis scales in the 
two panels of Figure |[ in reality, the weight of the peaks 
in the right hand panel is about 5 times higher than those 
in the left-hand panel) . These peaks form an almost dis- 
persionless band near £;=0, which then splits into two 
bands. These cross the Fermi energy separately and "re- 
combine" at k=n. The picture is somewhat unclear in 
that there is also a dispersionless "band" of low intensity 
peaks in the electron addition spectrum, which skims just 
above the chemical potential. We have scanned a vari- 
ety of hole numbers Nh and chain lengths N and always 
found either one of the two types of spectra in Figure 
|; more precisely, for N h /N=2/U, 2/16, and 2/18 we 
found a spectral function as in the left panel of Figure 
whereas for N h /N=2/12, 4/14, 4/16 and 6/16 we found 
that the spectral function looks like the right hand panel 
of Figure [|. This indicates that the very different shape 
of the two spectra does not originate from some spurious 
commensurate odering, which would occur only for one 
special value of N^/N; neither can it be due to a low 
spin-high spin transition, because all ground states un- 
der consideration are spin singlets. Rather, the only rel- 
evant quantity determining which spectrum is observed 
is the hole density, with the critical density S c for the 
crossover being l/7<<5 c <l/6. We thus have a concen- 
tration dependent phase transition between two ground 
states of very different nature. While details of the phase 
diagram will be reported elsewhere, we also note that 
for both high and low doping region the phase transi- 
tion to the 'anomalous' phases occurs for \t'\ « J (the 
precise value depends somewhat on J, being smaller for 
large J; both transitions can be clearly identified by a 
pronounced change of the electron momentum distribu- 
tion). The close relationship of the 'critical' t' with J 
also shows that the transitions are not driven by a defor- 
mation of the noninteracting band structure. We have 
observed the transition to the 'hole pocket' phase also in 
the Hubbard model with next-nearest neighbor hopping 
as U increases; this also will be reported elsewhere. 
We still investigate in more detail the properties of the 
high doping phase. Figure |^ shows the development of 
the spectral function for increasing hole concentration. 
A(k, ui) behaves very similar to that of noninteracting 
electrons, with the chemical potential progressively cut- 
ting more and more into an almost rigid "band" . The 
only exception is the dispersionless low intensity peaks 
in the addition spectrum, which seems "pinned" to the 
chemical potential. However, there is a very significant 
difference as compared to ordinary electrons. Inspection 
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of the sequence A//j=4, 6, 8 shows that removing two elec- 
trons from the system shifts the Fermi momentum by 
7r/8. On the other hand, to shift kp by this amount for 
spin-1/2 particles would require to remove 4 electrons 
(namely one electron/spin direction at ±kp). In fact, 
counting those momenta where PES and IPES spectrum 
have a strong low energy peak (such as 67r/8 in the 4 
hole case) as "half-occupied" , the number of momenta in 
the unoccupied part of the band equals Nh for Nh > 4. 
We thus arrive at the conclusion that the doping depen- 
dence of the spectral function in the high doping phase 
is consistent with holes behaving as spinless Fermions. 
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FIG. 5. Full single particle spectral function for a 16-site 
ring. Parameter values as in Figure ^j. 

In summary, we have studied the ID t—J model with 
an additional next nearest neighbor hopping integral t' . 
As a surprising result, the physics of the model seems to 
change completely already for moderate values of t' . In 
particular the data suggest a phase transition between 
phases of different Fermi surface volume. Thereby the 
low doping phase shows quite some similarity with the 
case of 2D, but also with doped t—J ladders: the "holon" 
excitation is overdamped, the single particle spectral 
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function j7j and the density correlation function [[lO 
consist of sharp low energy peaks and structureless con- 
tinua at high energies. Upon doping the spectral func- 
tion shows rigid-band behaviour [[L2| , with the chemi- 
cal potential cutting into the quasiparticle band seen at 
half-filling. The Fermi surface takes the form of a hole- 
pocket The location of the pockets is shifted to 



k=Q. As the density of holes exceeds a critical value 
S c w 0.15, there occurs a phase transition to a different 
ground state. The spectral function now shows a well de- 
fined band of sharp dispersive peaks, the doping depen- 
dence is consistent with the doped holes being spinless 
Fermions which gradually occupy this band. The rela- 
tion between electron density p e and Fermi momentum 
kp, being kp = |(1 — p e ) in the low doping phase and 
kp = 7r(l — p e ) for high doping, is therefore never con- 
sistent with the value for the Luttinger liquid, where it 
is k F = ^p e (it has to be kept in mind that the exact- 
diagonalization technique is hampered by the coarseness 
of the available momentum and energy resolution; how- 
ever, the Fermi momenta are quite unambiguously dis- 
tinguishable in the data and such qualitative results are 
not likely to be prone to finite-size effects). This result 
indicates a very profound reconstruction of the electronic 
structure as compared to the conventional Luttinger liq- 
uid. It also implies that these phases are not accessible 
(and hence have not been found previously) by standard 
methods for discussing ID systems, such as Bosoniza- 
tion |s| or renormalization group calculations Q. Both 
methods start out from the noninteracting Fermi points 
at ±§p e and attempt to construct effective Hamiltoni- 
ans for the low energy excitations around these. While 
the resulting ground state is qualitatively quite different 
from the noninteracting Fermi sea, it still "inherits" the 
period of its long range phase coherence, i.e. the corre- 
lation function (c\ a Cj_ a ) oc sm(k F \i — Anticipating 
that the Fermi momenta of the two phases of the t — t' ' — J 
model are indeed as discussed above, it is then the period 
of the long range oscillations of this correlation function 
are different, although the exponents for its decay are 
probably again consistent with "modified Luttinger liq- 
uids". 
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